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Abstract 

The issue of confinement and bose condensation is studied for gauge models 
of high-Tc superconductors. First the Abehan-Higgs model in (2+l)D, i.e., 
XY-model coupled to lattice gauge field with coupling g, is studied taking 
into account both the instantons and vortices. This model corresponds to 
integer filling of the bosons, and can be mapped to a dual superconductor. 
Our main resut is that the instantons introduce a term which couples linearly 
to the dual superconductor order parameter, and tend to pin its phase. As a 
result the vortex condensation always occurs due to the instantons, and the 
Meissner effect for the gauge field is absent, although is massive. This 
state is essentially the same as the confining phase of the pure gauge model. 
Away from integer filling, a "magnetic field" (the chemical potential of the 
bosons) is applied to this dual superconductor. Then the Higgs phase revives 
in the case of weak g and large x, where vortices do not condense in spite 
of the instantons. In the opposite case, i.e., strong g and small x, phase 
separation occurs, forming either microscopic patches or macroscopic stripe 
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domains of the Mott insulating state. 
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I. INTRODUCTION 



It has been established that the strong Coulomb repulsion between electrons is the key 
issue in the physics of high-Tc cuprates. Anderson proposed that this strong correlation 
gives rise to resonating valence bond (RVB) state, where the spin and charge are carried by 
spinous and holons, respectively |[T| . This phenomenon is called spin-charge separation, and 
has been subject to intensive studies. One way of formulating the strong electron repulsion 
is to exclude the double occupancy of the electrons on each site, and study the effective 
Hamiltonian within that restricted Hilbert space. Slave boson method is a useful tool to 
implement this constraint and fits the idea of spin-charge separation, where two species 
of particles, i.e., spinous (fermions f-^, fia) and holons (bosons 6|, bi), are introduced to 
represent electron operator Cj^, Cia- as f^-Q 

a. = hA- (1) 

Here an electron is represented as the composite particle of spinon and holon. At the mean 
field level, these two species of particles are supposed to be independent of each other. The 
mean field phase diagram is determined by two phase transition lines, i.e, the spinon pairing 
transition and holon condensation characterized by the order parameters A =< fiafja > and 
B =< bi >, respectively. In strange metal state both A and B are zero, while only A is 
nonzero in the underdoped "spin gap state" and only B is nonzero in the overdoped "Fermi 
liquid state" . The superconductivity is realized only when both A and B are nonzero, and 
the onset of the superconductivity is identified as the holon condensatoin in the underdoped 
region. However this simple-minded picture of spin-charge separation needs to be critically 
studied because the constraint is replaced by the average one in the mean field theory and 
the more appropriate treatment of this constraint might change the whole picture. This 
constraint can be taken care of by gauge field which corresponds to the fluctuation of the 
phase of RVB order parameter and Lagrange multiplier to impose the constraint ||5|-0. 
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Therefore the effective model is that of the spinons and holons coupled to the gauge field. 
It should be noted here that the gauge field represent the constraint and does not have its 
own dynamics. The effective action and the inverse of the coupling constant 1/g for the 
gauge field is generated only after integrating over spinons and holons. Therefore it is a 
highly nontrivial and crucial issue if the weak coupling perturbative analysis with respect 
to g makes sense or not. 

This issue is closely related to the confinement /deconfinement of the gauge field. |]5|,||,|^ 
The original model is defined on a lattice and the gauge field is compact, and it is well-known 



that the gauge field is confining in the strong coupling limit, i.e., large g. |[T0| A simplified and 
rough picture of the confinement follows. On the lattice one can define the gauge fiux h{p) 
penetrating each plaquette and the action is periodic with respect to h{j)) with period 27r. 
The simplest potential energy for h{p) is —g^^ cosb{p), and the kinetic energy of b{p) is given 
by where e is the electric field canonical conjugate to b{p). A phase transition is possible 
between two states. One is the extended "Bloch wave state" of b{p) where the tunneling 
events between different minima of the poriodic potential are driven by the large kinetic 
energy, i.e., large g. The conjugate field e is localized on the other hand, and the string 
of the electric field is formed when positive and negative gauge charges are inserted with a 
separtion R. This costs an energy proportional to R because of the finite string tension of 
the electric field. This phenomenon is called confinement. For small g, on the other hand, 
the periodic potential is large and b{p) is confined within one minima. One can replace 
the poriodic potential by the quadratic one ^b'^, which corresponds to the usual Maxwell 
Lagrangian. The Coulomb law is reproduced in this case, and the periodicity is irrelevant. 
The gauge field is deconfining in this case. In (2+l)D these tunneling events are represented 
by the point singularities of the gauge field configuration called magnetic monopoles or 



instantons. |Tl|] In the pure gauge model, the interaction between the magnetic minopoles is 
1/r, i.e., the Coulomb gas. Because the Coulomb gas in (2+l)D is in the screening phase, 
the gauge field is always confining due to these monoploes/instantons. 

When coupled to the holons and spinons, however, the gauge field becomes dissipative 



when these particles are integrated over, and the deconfining phase becomes possible in 



the strange metal normal state. |]T2[ However , in the presence of gapless fermion or boson 
excitations, the integration over them is in general not justifiable. The fermion part is 
better controlled because of the presence of a large energy scale, the Fermi energy, and 
techniques such as large N can be used to control the expansion. The boson part is much 
more problematic because bosons tend to condense in the bottom of their band and we are 
faced with a strong coupling problem of bosons and gauge fields. 

It has been argued that the strong inelastic scattering due to gauge fiuctuation suppresses 
the coherency and hence the ordering temperature, but the effects of the quantum fluctuation 
in the strong coupling limit is a difficult problem which remains unsolved. In this paper 
we hope that a duality mapping of a simplified version of the boson gauge field problem 
can shed some light on the issue. Summarizing the above, there are two crucial and related 
issues in gauge models of high-Tc superconductors, i.e., the Bose condensation and the 
confinement/deconfinement of the gauge field. 

In order to clarify these two pictures, we study in this paper a simplified model, i.e., a 
Higgs model coupled to U(l) gauge field defined on a (2+l)D lattice, which is an important 
model of broad interests both in condensed matter physics and high energy physics. fTBIJTl 
The action is given as 

S = -kY^ cos[A^6'(i) - ga^(z)] - - ^ cos[A^a^(i) - A^a^(z)] (2) 

link 3 plaquette 

where i is the lattice point and fi, v specify the direction in the (2+l)D lattice. The difference 
operator is defined as A^f{i) = f{i + jj) — f{i). Both 0{i) and a^{i) are compact and 
defined in the interval [0,27r]. q is the (integer) charge of the Higgs bosons and q = e 
(fundamental) for the holons while q = 2e for the spinon pairing in the U(l) gauge model of 
high-Tc cuprates 0]. Hereafter we take the unit e = 1 except for eq.(27) in section IV. The 
amplitude of the bose field has been fixed, but the vortex excitations are allowed due to the 
lattice and the compactness. Note that, in contrast to the high Tc problem, the gauge field 
dynamics is Maxwellian in the continuum limit. This model has been studied extensively. 



and the essential features are as follows. 



(1) In the limit g —>■ 0, the gauge fields are frozen out and we recover the XY- model. For 
(i + 1 > 2, we expect a phase transition at k = Kc- The ordered phase in d + 1 > 3 is 
characterized by an order parameter < e*^ >7^ 0. For g nonzero, however, this object 



is not gauge invariant and can no longer serve as an order parameter. [^] Instead we 
may fix the gauge to be, for instance, the unitary gauge ^ = and consider small 
gauge fluctuations. Then we have the Higgs mechanism where a term pga"^ appears in 
the action ( p^: the superfluidity density). The gauge flux correlation has the form 

< £ - ^1^) . (3) 

(2) The limit k = yields the compact QED model in 2 + 1 dimensions, which is known 



to be confining due to the appearance of instantons. |]TT| The instantons are singular 
configurations of a which act as magnetic monopoles, i.e., sources and sinks of magnetic 
fields. Writing a = + ainst as the sum of nonsingular and singular configurations, 
we have 

^inst = V X Oinst (4) 

For a given instanton density pM, we have the analog of Poisson's equation, 

V • 6inst = V ■ V X ttinst = 47rpM . (5) 

If the instantons form a gas, the density- density correlation function is given by 

< PM{k)pM{-k) >= -j^^—gi (6) 

in analogy with the more familiar density- density correlation function of a Coulomb 
gas, where Mq plays the role of the inverse screening length. Combining Eq. (4) and 
Eq. (5) we find 
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When combined with the standard transverse correlation from the nonsingular part 
Oo, we have 

, , r k^ky k^ki, Mq 

< 0^0. > - - + ^2 M^ + k^ 

k k 

Equation (7) shows that the electromagnetic field acquires a mass due to the gas of 
instantons. iTTH 



Thus we see that in both the g ^ limit and the k = limits, the gauge field is massive. 



This inspired Fradkin and Shenker [jT3[ to consider whether the two limits are smoothly 



connected to each other. Their conclusion is as follows: 

(1) When the charge is fundamental [q = 1), the strong coupling expansion argument shows 

that the Higgs phase with large k and small g continues smoothly to the confinement 
phase with large g and small k. 

(2) When the charge is not fundamental, e.g., q = 2, the Higgs and the confinement are 

different phases and can be distinguished by the forces between the test charges go = 
±1. 

In the bulk of this paper we will focus on the q = 1 case. The phase diagram for (3+l)D 
and q = 1 case has been determined at least qualitatively, and consists of two phases, i.e., the 
Higgs-confinement and Coulomb phases as shown in Fig. 1(a). In the Coulomb phase, the 
gauge field is deconfining and massless and the bose field remains disordered. However the 
phase diagram for (2+l)D and q = 1, which is the most relevant case to high-Tc cuprates, 
remains controversial [p!^p!3 |. One is that the XY-transition becomes first order once the 



coupling g to the gauge field is turned on, due to the Weinberg- Coleman mechanism [|16 . 
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This first order transition line terminates at some critical point. This resembles the vapor- 



liquid phase diagram. The other one is that a finite region of Coulomb phase exists. 



Since instantons dominate the physics for k = 0, the key question is whether instantons 
play an important role for general k and g. Vortex lines in the bose liquid carry unit 
fiux quanua and these can originate and terminate at instantons and anti-instantons. This 
is illustrated in Fig. 2. Einhorn and Savit discussed the free energy of the finite vortex 
segments, and concluded that the XY-transition remains at least for weak g. |]14| This leads 
to a picture of a finite region of the Coulomb phase mentioned earlier. Apart from the phase 
diagram, another issue is the behavior of the magnetic field correlation function in the region 



commonly labelled as confinement-Higgs after the work of Fradkin and Shenker. [T^ The 
point is that even though the gauge field is massive, the correlation function in the Higgs 
phase and the confinement phase are very different, being given by Eq. (2) and Eq. (7), 
respectively. The question is then whether the correlation function is closer to Eq. (2) or 
Eq. (7) in the confinement-Higgs region. As far as we know, this question has not been 
addressed. In this paper we include the effect of instantons for general n and g but we come 



to quite a different conclusion than Einhorn and Savit. 11^ The phase diagram we propose 
is shown in Fig. 1(b). The line = is an isolated line with an isolated XY-transition. We 
also calculate the magnetic field correlation function and show that it takes the confinement 
form of Eq. (7). Thus we conclude that the Higgs-confinement phase is better described as 
the confinement phase. This is the main result of this paper which we discuss in Section II. 

It is known that the model given by Eq. (1) describes a Bose gas on a lattice where the 
density is an integer. Thus the confinement phase may be understood as a Bose Mott insu- 
lator. In Section III we extend the discussion to the incommensurate case, when the density 
is not an integer. We find two possibilities. The instantons may become irrelevant and the 
system becomes a superfiuid (Higgs phase), or the system may break up into domains. In 
Section IV we briefiy address the case q = 2 and also the case of two kinds of bosons with 
opposite gauge charges, which arises in the SU(2) formulation of the t-J model. 
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II. THE ABELIAN-HIGGS MODEL 



We start with the continuum action for eq.(l 

1 



S 



d X 



1 



-k(V^(x) - a{x) - A{x))^ + —(V X df 
/ Zg 



(9) 



Here we have introduced the external electromangetic field A^, which is put to be zero for 
the moment. We allow the singular configurations of 6 and a corresponding to the vortex 
and instanton, respectively. As for the vortex, 6{x) is devided into the single- valued part 6q 
and the multi- valued vortex part By d.s 6 = 6q + 6v, and V x VOv/^tt = jy is the vortex 
current. As for the instantons, a = ao + flinst and V • V x ainst/47r = pM is the instanton 
density as defined in Eqs. (4) and (5). We next derive a continuum version of the duality 
representation of this problem. The duality representation is a powerful tool which allows 
us to discuss the strong coupling limit ( large g ) and is particularly useful for the present 
problem. Introducing the Stratonovich-Hubbard field J representing the boson current, the 
action becomes 



S 



d X 



—.P + ij- iweix) - dix)) + —(V X df 

2k 2g 



(10) 



and after integrating over Oq, we obtain V ■ J = corresponding to the conservation of 
the boson current. To enforce this conservation law, we introduce the vector potential c to 
represent J as J = V x c. Then after partial integration the action becomes 



S 



d^x 



^(V X C-)2 + IC- {27CJy - 6(x)) + ^{bf 

2k 2g 



(11) 



with 6 = V X a. After integrating over b, we obtain 



S 



d-^x 
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— {Vxc)' + ^{c)' + t27rc-jv 



(12) 



.2k' ' 2 

This is essentially equivalent to the action obtained by Einhorn-Savit in terms of the lattice 



formulation It is noted that if we view c as a gauge field coupled to the vortex current, 
the gauge symmetry is broken in Eq.(ll), which corresponds to the non-conservation of 
the vortex current, i.e., V ■ jv can be nonzero, due to the instantons. After integrating 
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over the field c in Eq.(ll), the partition function Z is given by the integral over the vortex 
configurations {jy} including both the vortex loop and open vortex segments the two ends 
of which are instanton and anti-instanton . 

^ = ^0 XI expE -4vr^/«jy^(j)L)^^(j - k, rr?)3vv{k)\ (13) 
Civ} 

where Df^i^{j — k, m?) is the propagator of the field c with the mass m? = Kg. This propagator 
is given as 

A„A, 



where D{j — k,m?) satisfies 



^fl'—^U 



D{j-k,m^) (14) 
j 



{-Al + m')D{j-k,m') = 5,k. (15) 

The propagator decays exponentially as D(j — k, m^) ~ and the interaction between 

the vortex segments are short range. 

Let us start with a qualitative estimate of the free energy of the vortex loop/segment 
regarding the partition function as that of a classical statistical mechanics. We first repeat 



the argument of Einhorn and Savit. Let L be the length of the vortex loop/segement. 



Then the free energy is the sum of the energy cost and the entropy as fT^ 

Fioop = 4'K\D{0;m^)L - ln/2^ 
i^segment = ^7r^KD{0; m^)L - ln(/i^L^) + 25i„st (16) 

where S'inst = 47r^(7~^D(0; m^) is the action for a instanton, and /i is a number of order 
unity which depends on lattice and dimension. The main difference between closed loop 
and open segment is that the two instanton action is added and entropy is enhanced by the 



factor U' (7 > 0) in the case of open segment |]T^. However the leading L-linear terms are 
the same for both of closed loop and open segment, and Einhorn-Savit concluded that the 
proliferation of the vortices, i.e., the appearance of the infinite length loop/segment, occurs 
at 47r^fi;D(0; m^) = Injl. This is the XY-like phase transition viewed in the duality picture. 
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However the above consideration neglects the possibihty that the long vortex loop/segment 
are cut by the instantons into small pieces, as shown in Fig. 2c. Let us consider that the 
total length L is cut into n pieces of open segments. The free energy in this case is 

n) = (ATT^KDiO; n?) - fx)L - -fn ln(L/n) + 2Sinst"- (17) 

Minimizing this with respect to n, we obtain the length of the pieces as L/n = 6^+^'^'°=*/''', 
which is finite as L oo. Then due to the finite density instantons, the infinite length loop 
or segment does not appear. Instead the vacuum is full of finite size open segments and 
closed loops of vortex. The total length of these loops/segments are infinite as the sample 
size L ^ oo when 4it'^k,D{0; m?) — /i < 0, but they are all finite size. Therefore we conclude 
that the phase diagram of (2+l)D Abelian Higgs model is given by Fig. 1(b), i.e., all the 
interior constitutes a single phase. Analogy of the phase diagram with the Ising model is 
useful here. In the duality picture, i.e., regarding the vortex field as the order parameter, k 
can be regarded as the temperature T, and the instanton density e"™"^*/^ is the magnetic 
field H. Then it is natural that the ordered state aXl/ g = oo, i.e., if = 0, is the isolated line 
while all the other phase diagram is connected to the high temperature symmetric phase. 
Actually this analogy becomes more clear when one consider the path integral formulation 
of the the Ising model under magnetic field H [l^ll^ • The partition function is represented 



by the sum over the closed loop and the open segment ended at the magnetic field vertex 
H. This is exactly similar to the present case where the vortex segment terminates at 
instantons except that the vortex loop/segment has a direction and instanton has the ± 
topological charge. Now what is the nature of this single phase ? One crutial question is if 
the Meissner effect for a remains or not. It should be noted that the magnetic field h is tied 
to the vortices. In the case of closed loop, however, the net magnetic field h is zero. The 
open segment, on the other hand, is a magnetic dipole which has net h. Therefore once the 
instanton fugacity is nonzero and there are vortex open segments, the magnetic field h can 
penetrate into the sample, i.e., the Meissner effect disappears and the gauge field h becomes 
massless. This corresponds to the "dielectrics" of the magnetic charge and is illustrated in 
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Fig. 3b. The instantons and anti-instantons appear to be bound into pairs. However this 
state is not stable because once the Meissner effect is gone, the confinement between the 
instanton and anti-instanton also disappears, and the magnetic dipole is liberated into free 
magnetic charges. This is the "metal" of the magnetic charge, and again the magnetic field 
b can not penetrate into the sample due to the screening. This is shown in Fig. 3c. This 
massive gauge field b corresponds to the confining phase of the pure gauge model discussed 
by Polyakov Therefore the state in the interior of the phase diagram in Fig. 1(b) 

continues smoothly not to the Higgs phase but to the pure gauge model, and should be 
called "confinement". 

In order to substantiate the above consideration, we now go to the second quantization 
formulation of the vortex system. Let us first devide the gauge field c into the transverse 
and longitudinal parts, i.e., c = c± + V0. Then the action eq.(ll) is written as 



S= (fx 



^(V X cl)' + |(cx)' + |(V0)2 + i27rcx ■Jv + i<p-pM 



(18) 



where we have performed a partial integration and identified V- jy with the instanton density 
Pm- We can view Eq. (18) as describing world-lines of vortex-particles which are coupled 
to the gauge field by the i2TTc± ■ jy term. The vortices may be created and annihilated at 
instantons located at Xi ■ ■ ■ x„ and anti-instantons located at ?/i ■ ■ ■ y^. Alternatively, we can 
write the action in terms of the second quantized vortex field ipv The action is given by 

Sv= f d'xl^iV X cl)2 + |(c^)2 + |(V0)2 + ^lll-KiV + tc^f + M^]^v + u{^Wv) 
J iZK 2 2 2 

(19) 

where = 47r^/tZ}(0; m^) — Infl, and u represents the short range repulsion between the 
vortex segments. The gradient term comes from the extra cost of the action when the 
vortex line deviates from the straight line. This step is standard in the duality mapping. 
The novel feature of creation and annihilation of vortices can be included by summing over 
all instanton configurations as follows, 

dxi ■ -dxn f dyi ■ -dy,^ 



Z = J Di;iDi,yDc^D<P J ^'^j / 



n,m=0 ' 
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where z is the fugacity of the instantons which is roughly given as 2; ~ e~'^'"=*. n,m are 
the number of instantons and anti-instantons, but only the term n — m survives when one 
integrates over ipv,'4^v- -^^ the usual Coulomb gas mapping, the summation over n,m in 
eq.(20) can be done, and our final result for the action when recovering the original gauge 
field a is given as 



(fx 



7^(V X c±f + ^{V X af -ic-Vxa 
2k 2g 



(21) 



where it should be noted again that c = cj_+V0 and a = ao+Oinst- The gauge transformation 
ipv ^v^^'^j V'y ~^ ^ve~^'^ eliminates the exponential factor in the 2;-term, and also replaces 
cj_ by c in the minimal coupling term. We have 

J DijjlDijjvDcDae-^ (22a) 
S ^ [ d^xl ^{V X c)^ + ^(V X df - ic-V xd 

V) 

-z(4 + V'y)[ (22b) 



Except for the last term, Eq.(22b) is the standard duahty representation of the abehan 
Higgs model. It is useful to recall that vortices in the field z/'y corresond to world fines of the 
original bosons. Condensation of the vortex field i\}y means the absence of Bose condensation 
and vice versa. The last term in Eq. (22b) represents the effect of the instantons and is 
the main new result of this paper. We note that it takes the form of an external field 
coupled to the vortex field V'y. This is analagous to the Josephson coupling to an external 
superfiuid with an order parameter z. The external order parameter will induce a nonzero 
order parameter < ■0V >^ 0, even when > 0. [Note we have fixed the gauge and the 
nonzero order parameter is in the particular fixed gauge.] In the first quantized picture of 
vortex loops and segments, this can also be understood as follows. Consider the correlation 
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function C{x,x') =< tpy{x)tpv{x') >■ In the absence of instantons, the points x and x' are 
connected by a vortex hne and long-range order is possible only when an infinite vortex line 
has zero energy, i.e., < 0. However, with finite z, an instanton and an anti-instanton 
appear near x' and x and create two finite segments, so that C{x,x') reaches a finite value 
even as the separation between x and x' goes to infinity. 

Recall that in the duality picture, < ipy >7^ means that the original boson is not 
bose-condensed. Thus we expect that the effect of the z^ipy + V'y) term is to destroy the 
Meissner effect of the original Abelian-Higgs theory. We check this by an explicit calculation 
of the gauge field correlation function. Let us represent ipv as ipv = V'oc*'^. Then the action 
for the vortex field becomes 

S^orte. = J d'^X^ [Kll^oC^(p + C}^- ^^^0 COS (^] (23) 

which is the sine-Gordon model in (2-|-l)D. According to the analysis of the pure gauge 



model by Polyakov [|Tl|], the fugacity z is always relevant and (/9-field is massive. Replacing 
the cos-term by the effective quadratic term as 

^vortex = J d^x^[K^l{V^ + c)^ + 2z,s^o¥>^] (24) 



which gives the mass mo = \j2zes/ {Kipo) of the y^-field corresponding to the screening. After 
integrating over we obtain the effective action for the gauge fields c, a and the external 
electromagnetic field A as 

+ ^((^/.z. - - i^i^Ai^M) + A^{q))cy{~q) + c^(-g)(a^(g) + A^q))] 

After integrating over c, we obtain the propagator of the gauge flux 6 = V x a as 



(25) 



1 
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It should be noted that while the pole disappeared and the gauge field is massive, there 
is no Meissner eff'ect because Eq. (26) is not proportional to as in Eq. (2). Rather, the 
small q behavior is essentialy the same as that in the confining phase in the pure gauge 
model shown in Eq. (7) ||lT|. After integrating over the h field, the effective action for the 
external e.m. field is ~ ( V x A)^, which means that the system is insulating. This is perhaps 
not surprising if viewed in the strong coupling limit g oo. Then the gauge field does not 
have its own dynamics and serves to impose the constraint of integer occupation at each 
lattice site. The bosons are just frozen into place on each site, resulting in an insulator. 
This Mott insulator phase appears to extend to include the entire phase diagram, as shown 
in Fig. 1(b), with the exception of the line ^ = 0. 

At finite temperature T, the imaginary time axis becomes finite, i.e., [0,(3 = 1/T], in 
eq.(23). Therefore eq.(23) desciribes the sine-Gordon model in 2D in the long wavelength 
limit. Therefore we expect the KT transition, i.e., confinement-deconfinement transition, at 
some critical temperature T^. 



III. BOSONS WITH NON-INTEGER DENSITY 

The Abelian Higgs model in eq.(l) corresponds to the case of integer boson filling for 
each site. The deviation from it is taken care of by introducing the winding number of the 



phase 6 along the time direction |T9|. In the dual picture, the boson density is represented 
by the z-component of the gauge fiux V x c, and the deviation from the integer filling is 
represented by adding the term — yu(V x cj^ to the action eq.(22). Here the chemical potential 
fi acts as the magnetic field. Therefore the chemical potential and the vortex condensation 
compete with each other as in the case of the superconductor in a magnetic field. The 
only new aspect here is the instanton term z{tpv + i^v)- When fi = 0, this term induces 
the vortex condensation < ipv > even if is positive and large. For /i 7^ 0, we may be 
tempted to consider an Abrikosov vortex state (of the vortex field ipv) in analogy with type 
II superconductors. An ordered array of such vortices correspond to a Wigner crystal of 
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boson. [19] However the phase of ipv changes by 27r around each vortex and we cannot gain 
the Josephson energy from the term zlipy + i^v)- We conclude that the Wigner crystal is 
suppressed by the instantons. A second idea is to consider the analogy of the intermediae 
state in type I superconductors, where the stable configuration is the laminar structure 
. In type I superconductors the surface energy between the normal and superconducting 



regions is positive. The surface energy is proportional to the size along the 2;-axis {(3 in 
the present context) and the spacing of the laminar structure is macroscopic in size. In the 
appendix, we perform a Ginzburg- Landau calculation of the surface energy, and find that 
in contrast to usual type I superonductors, it is negative in the case > 0, i.e., when the 
superconductivity is induced by Josephson coupling. This implies that the straight interface 
is unstable, and the laminar phase will break up. One possibility is that the system breaks 
up into patches where < ipv 0^ separated by regions where < ipv >= 0. (This can be 
viewed as the complement of the Abrikosov vortex state.) The order parameter can be real 
in each patch, gaining an extensive Josephson energy from the term —z^ipv + i^v)- The 
magnetic field V x c can penetrate the normal region and partially penetrate the patches. 
This state can maximize the surface energy gain for a fixed patch area and the patches will 
form some ordered structure. In the original boson representation, the absence of long-range 
order in ipv means that the bosons form a superfiuid with Meissner effect. The instantons 
become irrelevant in this sense, in contrast to the /i = case. The effect of the instanton 
is to cause a periodic modulation of the boson density (corresponding to the modulation 
of the magnetic field V x c by the patches). This modulation is weak for k large (M^ 
positive and large) and grows with decreasing k. An ordered array of patches lead to a kind 
of incommensurate order. For very small k, < and there is a strong tendency for 
the order parameter ipy to form in the dual picture. In this case the interface energy may 
become positive and we cannot rule out a laminar picture. In the original boson picture this 
corresponds to stripes of superfiuids separated by Mott insulators. The transition between 
the stripe phase depends on details of the parameters and we have not attempted to work 
it out quantitatively. However, since stripes occur only for positive interface energy in our 
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model, the stripe size in expected to be macroscpic , by ananogy with the laminar phase of 
type I superconductors. Microscopic stripes might occur when < and other interactions 
such as long range Coulomb forces and commensur ability energy are introduced. We have 
not examined this issue in this paper, and this is left for future studies. 



Up to now, we consider the charge e bosons coupled to the gauge field. However other 
types of order parameter appear in gauge models of high-Tc superconductors, which is the 
subject of this section. 

First we consider the case of g = 2e, which corresponds to the spinon pairing order 
parameter coupled to the U(l) gauge field. Here we take the unit where 2e — 1 and the flux 
quantization is reduces to half. Then the instanton becomes the end point of two vortices, 
which modifies the ^-term in eq. (22b) as 



It is noted that 2;-term is the quadratic term and does not necessarily enforce the condensa- 
tion of -01/, Therefore two possibihties arises in this case. 

(I) Single vortex condensation, i.e., < ipy >7^ 0. In this case the quantized charge, i.e., 

the integral of (V x c±)z is 2e and the single charge e can not appear. Therefore the 
charge e is confined. 

(II) Vortex pair condensation, i.e., < '4>v'^v while < V'y >= 0. In this case the 
quantized charge is reduced to half, i.e., e. Therfore the confinement of the charge e 
does not occur. 



IV. OTHER ORDER PARAMETERS 




(27) 



17 



Fig. 4 shows the phase diagram for q = 2e, where the above two possibihties correspond 
to I and II, respectively. As for the small fluctuation of a is concerned, there occurs no 
Meissner effect in both phases, although a is massive due to the confinement. Therefore 
both phases is better called confining phase. What distinguished these two phases is the 



discrete Z2 symmetry. |TB| It has been discussed that the limit k = 00 is the Ising gauge 
model, which shows confinement and deconfinement transition at some critical value of 
9 = 9c 0. 

Secondly we study the two species of bosons bi, 62 coupled to the gauge field with opposite 
charges e and — e, respectively. This situation occurs in the staggered flux state of an SU(2) 
formulation for underdoped region. [21] In this case the dual model is given by 

S = f d'x{^[{V X el f + (V X c,Y] + 1(V X af 
J 2k 2g 

- 2cl ■ V X (a + Ai) -ic2-V X (-a + A2) - i/i(V x (cl + cl))^ 



z{4ii'2 + V^IV^i) + + {42i^2f] + 2«;(^iVi)(^|^2)} 



(28) 



where q {i = 1,2) is the gauge field representing the boson current of bi, ipi is the corre- 
sponding vortex field and Ai is the test field coupled to it. The real electromagnetic field 
corersponds to Ai = A2 = A. It is noted here that again the instanton term (2;-term) is 
quadratic in ip^s and cannot induce nonvanishing i/ji and ip2 when z is not large enough. The 
Bose condensation of bi and 62 should occur in this case.. When g and z is large, which is 
relevant to the high-Tc problem, the amphtudes of both ipi and 'ip2 are induced and we write 
'0j = V'oe*'^*- Here the singluar vortex configuration of (fi is allowed, which corresponds to 
the original boson. Then the effective action for the phase field (fi is given by 

■1 



5-, 



eff. 



-K^liVifi + ciY + {Vip2 + C2)'] - 2z4 cos(y^i - (P2) 



(29) 



Here we define the symmetric and antisymmetric parts as 



^1 = (fis- ^^a 



(30) 



and Csi Ca in a similar way. Then eq. (29) is wirtten as 



1 1 

-KlPl[2{Vips + Csf + -(V(^„ + Caf] - 2ziPlcOs{ipa) 



(31) 



where only the antisymmetric part is coupled to the instantons. Therefore the actoin for the 
antisymmetric part is the same as that in eq. (23), i.e., is fixed by the 2;-term, and the 
vortex of (/^a-field is forbidden. There is no Meissner effect for the gauge field a, although it 
is massive as in the pure gauge model. This corresponds to the binding or confinement of 
the two species of bosons hi and &2 because (V x Ca)z is the difference between the boson 
densities of hi and 62- This means that the single Bose condesation is suppressed. The 
boson pairing condensation, on the other hand, is not disturbed by the 2;-term. Therefore 
when the field is disordered, we have the boson pair condensation and finite superfulidity 

— * — * 

density ps- Then the effective action for the test fields Ai, A2 is given after integrting over 

c-fields as 



(Tx 



Ps{Ai + A^f + x„(V X {Ai - I2))' 



(32) 



where Xa is a diamagnetic susceptibility of the antisymmetric part. The the system show 
the Meissner effect only for the symmetric test field Ai + A2. Therefore the system shows 

— * 

the Meissner effect to the external electromagnetic field A. 



V. CONCLUSIONS 

In this paper, we studied the interplay between the confinement and the condensation 
of the order parameters. At integer-filling of the bosons with charge e, there is only one 
phase in (2-|-l)D with the XY-transition isolated only on the hne g — The nature of 
this so-called Higgs-confinement phase is the same as the confining phase of the pure gauge 
model, and no Meissner effect for the gauge field occurs. This is because the instantons act 
as ordering field for the vortex condensation. For non-integer filling, Bose condensation is 
recovered for weak coupling g. However, the Bose condensation and confinement compete 
with each other, and this competition leads to phase separation for strong coupling. 
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In this paper we have focused our attention to the problem of bosons coupled to gauge 
fluctuations. It is only a first step towards addressing the problem which arises out of the 
gauge theory formulation of the high Tc problem, which involves both fermions and bosons 
coupled to gauge fields. Nevertheless, we would like to put the present work in the context 
of the high problem and attempt to draw a few inferences. The effect of the fermion is 
two- fold. First, if the gauge field is confining, it allows the possibility of confining fermion- 
antifermion pairs to form spin excitations and confining fermion and boson to re-constitute 
the physical electron. The former is believed to happen in the half-filled case where the AF 
ordering may be described as chiral symmetry breaking and confinement of Dirac fermions. 



2^ , |25[| Secondly, the presence of massless Dirac particles changes the dynamics of the gauge 
field and, in general, it would not take the Maxwellian form assumed in this paper. We can 
divide the doping region of the phase diagram into three regimes: 

(i) Doping into the AF {x <^ 1). Here the starting point is a vr-fiux phase for the fermions 
with a Dirac spectrum. The gauge propagator is proportional to \/q^ instead of in 
the Maxwell theory. P The instantons have logarithmic interaction and undergo a 
Kosterlitz-Thouless transition in the 2 + 1 dimension as a function of A^, the number 
of fermion flavors. It is believed that the physical case of = 2 lies on the disordered 
side of this transition, so that the instanton gas behaves as free gas (as opposed to 
instanton anti-instanton bound pairs). Since our consideration is based on assuming 
the existence of the free instanton gas, this would be the case where our consideration 
has the best chance of being applicable. Nevertheless, we still have not included 
the possibility of bosons combining with fermions to form physical holes in an AF 
background. This would correspond to the formation of small Fermi liquid pockets 
in a reduced Brillouin zone. Leaving this possibility aside, we can conclude from 
the results of Section III that instantons suppress the formation of a Wigner crystal 
of doped holes. Furthermore, the possibility of phase separation into microscopic 
patches is interesting, in that it suggests incommensurate structures which appear 
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experimentally in this part of the phase diagram. However, the superfluid state that 
appear in our picture does not appear to resemble d-wave pairing, as long as the 
fermions remain confined in the AF state. It is also interesting that phase separation 
into larger scale laminar domain is a possibility. Finally in the SU(2) formulation, the 
result of Section IV suggests the possibility of bosons forming a pairing state, leading 
to a co-existence of superconductivity and AF. 

(ii) Underdoped region. Here the normal state is the pseudogap state which is described as 
(i-wave pairing of fermions or a staggered flux phase. Again initially the fermion 
spectrum is Dirac and the gauge propagation is proportional to and it is not 
clear that the present paper is applicable. Nevertheless, we can ask whether the 
low temperature phase is a conflnement phase where instantons are free and play an 
important role. There are three possible scenarios for the onset of the low temperature 
superconducting phase. The flrst is a binding of fermions with bosons to form physical 
quasiparticles. Since the fermions are already paired, a superconducting state appears. 
This possibility is clearly beyond the scope of the present work. The second and third 
possibilities are the Bose condensation of single bosons in the U(l) formulation, or the 
pairing of two kinds of bosons in the SU(2) formulation. |^ The latter problem is 
treated in Section IV. What we learn from the present study is that instantons tend 
to suppress Bose condensation when the coupling constant g is larger. Furthermore, 
instantons favor the binding of the two species of SU(2) bosons to form pairs which 
then condense, leading to a d-wave superconductor ground state. 

(iii) The overdoped region. Here the high temperature phase is the strange metal phase 
and it has been argued that it is a de-conflning phase due to dissipation in the gauge 



field dynamics. |12| The low temperature Fermi liquid phase is best described as a 
confinement of fermions and bosons. These are clearly outside of the scope of the 
present paper. 
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APPENDIX A: SURFACE TENSION 



In this appendix, we show the calculation of the surface tension between the normal and 



superconducting regions of the dual superconductor following ref. ||2^. The free energy of 
the dual superconductor is given by 

1 ^ .o 1, 



F = j dh [—(V X c^f + ^-[-/^(V + ic^f + M^]i;v - z{^v + 4) + ^i'^W' 



2 

V) 



(Al) 



where k'~^ = k,+ {dn/dfi)~^ ( dn/dfi: the charge compressibility ), and the classical ( time- 
independent ) configuration is assumed. The Ginzburg-Landau equations are obtained by 
taking the variation with respect to Sipy and 6c. 

K 

— {-iV-CA_fiJv + —i^v + u\iJv?i)v = z (A2) 
V X (V X cl) = ft'Jy (A3) 

where 

3v = ^(^y^VV^y - (V^v^) ■ iPv) - A\^v?c^ (A4) 

Now we consider the case of instanton driven dual superconductivity. Namely the 2;-term 
is the driving force of the vortex condensation and > . Then we assume is large 
enough and u{il)yil)vY term can be neglected. In the absence of the magnetic field c±_ , 
'ipv = ipvo = , and the free energy measured from that in the normal state F^o is given 
by 

F - F„o = -V^ = -V—, (A5) 
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where He is the thermodynammic critical field and V is the volume of the system. Assume 
that the interface between the superconducting and normal regions is localized near x = , 
and x > region is superconducting. Physical quantities depend only on x, and we choose 
the Coulomb gauge V • c_l = 0. Therefore dxC^^ — 0, and we put c^^ — . The boundary 
condition is 

dc I 

(V X cl)^ = = //c, 'i/'y = -oo 

(V X cl)^ = ^ =0, = -0^0 X ^ +00 (A6) 

Here we introduce normalized quantities. 

X — x/X, 

i> = i>v/i>vo: 

c^c^y/{H,X), (A7) 
where A is the penetration depth and is given by A = (AKk'Z'^/M^)'^/'^. The correlation 



length ^ is given by ^ = ^JK/\M'^\, and we define the ratio 77 = A/^. Using these normahzed 
quantities, the GL equations become 

^" = r/2[l - (c^ + c" = alj^, (A8) 

where ip" = (fip/dx'^ etc., and the boundary condition is 

c' = 1, ip' = a; — >• —00 

c' = 0, -0 = 1 +00 (A9) 

It can be easily shown from eg.(A8) that 

+ (c' + - 2^^ - c'2 = -1 (AlO) 

The surface tension q;„s is given as follows. First define the free energy density / under the 
magnetic field H as 
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/ = / 



HB 



K' 



(All) 



Then ans is given by 



/oo _ _ 

dx{f - /„) 
-oo 



dx 



2k' 



. "n/'i^.^i/ |2^,^^/ |2 // , lU HciVxc±), 2z'^ 
+ +c±|'0y| ) + ^IV'yl - ziijjv + i^v) -. h — 



K' 



dx[{c' - If + -{il^'f + + 1)V;2 - 2^], 



2k' J-oo " ■ r]^ 
where /„ is the / in the normal state. Using eg.(AlO), 

A/72 oo 



M2J 
(A12) 



^ ^ / dx(^{(^ - 1) 

K J —oo 



(A13) 



Because the normahzed magnetic flux density c' is < c' < 1 in the interface region, 
the integral in eq. (A13) is negative and < 0. Therefore we conclude that the surface 
tension is negative in the instanton-driven dual superconductor. 



24 



REFERENCES 

[1] P. W. Anderson, Science 235, 1196 (1987) ; The Theory of Superconductivity in High-Tc 
Cuprates Superconductors, ( Princeton University Press, 1997). 



[2] 
[3] 
[4] 
[5] 

[6] 
[7] 



[9 
[10 

[11 
[12: 

[13: 

[14 
[15 



G.Baskaran, Z. Zou, and P.W.Anderson, Solid State Commun. 63, 973 (1987). 
G. Kotliar and J. Liu, Phys. Rev. B38, 5142 (1988). 

Y. Suzumura, Y. Hasegawa and H. Fukuyama, J. Phys. Soc. Jpn. 57, 2768 (1988) 

G.Baskaran and P.W.Anderson, Phys. Rev. B37, 580 (1988); G.Baskaran, Phys. Scr. 
T27, 53 (1989). 

L. B. loffe and A. I. Larkin, Phys. Rev. B39, 8988 (1989). 

N. Nagaosa and P. A. Lee, Phys. Rev. Lett. 64, 2450 (1990); P. A. Lee and N. Nagaosa, 
Phys. Rev. B46, 5621 (1992). 

R. B. Laughhn, in Proceedings of the Inauguration Gonference of the Asia-Pacific Genter 
for Theoretical Physics, Seoul National University, Korea 4-10 June 1996, ed. by Y. M. 
Gho, J. B. Hong, and G. N. Yang (World Sci., Singapore, 1998) 



R. B. Laughlin, |cond-mat/9709195 



For a review on lattice gauge theory see e.g. J. B. Kogut, Rev. Mod. Phys. 51, 659 
(1979). 

A. M. Polyakov, Gauge Fields and Strings (Harwood Academic Pub., 1987). 

N.Nagaosa, Phys. Rev. Lett. 71, 4210 (1993). 

E. Fradkin and S. Shenker, Phys. Rev. D19, 3682 (1979). 

M. B. Einhorn and R. Savit, Phys. Rev. D17, 2583 (1978); ibidDlQ, 1198 (1979). 
S. Elitzur, Phys. Rev. D12, 3978 (1975). 



25 



[16] B. I. Halperin, T. C. Lubensky, and S. K. Ma, Phys. Rev. Lett.32, 292 (1974). 

[17] P-G. de Gennes, Scaling Concepts in Polymer Physics, (Cornell University Press, 1979) 
Chaps. X, XI. 

[18] J. des Cloiseaux, J. Phys. (Paris) Lett. 36L, 281 (1975). 
[19] M. P. A. Fisher and D. H. Lee, Phys. Rev. B39, 2756 (1989) 

[20] P-G. de Gennes, Superconductivity of Metals and Alloys, (Addison- Wesley, 1989) 
Chap.2. 

[21] X.-G. Wen and P. A. Lee Phys. Rev. Lett. 76, 503 (1996); 80, 2193 (1998). P. Lee, N. 
Nagaosa, T.K. Ng, and X.G. Wen, Phys. Rev. B57, 6003 (1998). 

[22] J. B. Marston, Phys. Rev. Lett.64, 1166 (1990); Phys. Rev. B42, 10804 (1990). 

[23] L Affleck, Z. Zou, T. Hsu, and P. W. Anderson, Phys. Rev. B38, 745 (1988). 

[24] E.Dagotto, E.Pradkin, and A.Moreo, Phys. Rev. B38, 2926 (1988); E. Dagotto, J. B. 
Kogut, and A. Kocic, Phys. Rev. Lett. 62, 1083 (1989); 

[25] D. H. Kim and P. A. Lee, Annals of Physics 272, 130 (1999), and references therein. 

[26] E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics - Part 2, (Pergamon Press, 1980) 
p 184. 



26 



Figure captions 



Fig. 1. Phase diagrams for Abelian Lattice Higgs model with fundamental charge in (3+l)D 
(a) and (2+l)D (b). In (b) the phase transition is isolated along the XY-model line, i.e., 
g — 0. In the confinement phase, the gauge field is massive due to instantons, althogh no 
Meissner effect occurs. 

Fig. 2. (a) Vortex loop, (b) a vortex segment, and (c) vortex segment broken into pieces by 
the instantons. A vortex segment terminates at the instanton and anti-instanton. 
Fig. 3. Three possible state for the magnetic charges, (a) The vacuum state where only 
vortex loops of finite size exist. The gauge flux b can not exist inside the sample, i.e., 
the Meissner effect occurs, and this is the superfluid state of the original bosons, (b) The 
dielectrics with finite size vortex segments and loops. The magnetic field b can penetrate into 
the sample, and the Meissner effect vanishes. The string tension of the vortex then becomes 
zero, and the confinement of the instanton and anti-instanton disappears. Therefore this 
state is unstable to (c). (c) The metallic state of the magentic charges. The metallic 
screening prevents the gauge field from penetrating into the sample. This is the confining 
state of the gauge field. 

Fig. 4. Phase diagram for q = 2e. In region I the single vortex condensation occurs and 
charge e is confined. In region II only vortex pairs condense and charge e is not confined. In 
the limit k — oo, the model is reduced to Ising gauge model, which shows a phase transition. 
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